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Abstract 

We consider diatomic systems in which the kinetic energy of the electrons is treated in a simple 
relativistic model. The Born-Oppenheimer approximation is assumed. We investigate questions 
of stability, deducing bounds on the number of electrons, the binding energy AE'f, and the 
equilibrium bond distance Rq. We use a known localization argument adopted to the present 
relativistic setting, with particular consideration of the critical point of stability, as well as the 
recently proved relativistic Scott correction. 



I. INTRODUCTION 



In this paper we study the stabihty of diatomic systems in the Born-Oppenheimer 
formulation. The kinetic energy of the electrons is treated in a model which accounts for 
relativistic effects. We employ a known localization argument applied by Solovej [l| to prove 
a lower bound on the number of electrons. This involves establishing an estimate of the 
equilibrium bond distance Rq as a function of which requires particular consideration 
in the relativistic case due to the critical point of stability for large values of the nuclear 
charges. Bounds on the binding energy AEb and on Rq are deduced using the recently 
proved relativistic Scott correction j2|. 

We introduce positive numbers Zi, Z2, Z = {Zi, Z2), which correspond in a physical picture 
to nuclear charges of a diatomic molecule, and denote the nuclear coordinates as Ri, R2 & 
M^. We consider furthermore a set of electrons and introduce the electron coordinates 
Xi G M.^, i = 1, . . . , N. The coordinate origin is placed at |(-Ri + i?2)- Let Ri = R/2, 
R2 = —R/2, where R = Ri — R2, furthermore R = \Ri — R2\. We are working in units 
where m = e = h = 1 (m: electron mass; e: electron charge; h: Planck's constant). In 
this unit system, the distance is measured in Bohr radii Oq = fi^/me^. The Hamiltonian 
employed in the present context is then 

H(N, Z, fl; a) ^ f (V-a-A...- - a-) - f . ^) 

i=l 1=1 ^ ' 

El Z^ Zo 

I _ I 

Kj ' 

= + Hen + Hee + Hnn- (l) 

a denotes the dimensionless fine-structure constant. It can be viewed as a quantifier of the 
relativistic correction and is in our unit system defined as a = c~^, where c denotes the speed 
of light. In the nonrelativistic limit, we have a — )■ 0. We assume the Born-Oppenheimer 
approximation, i.e., fixed nuclei positions -Ri, R2 corresponding to infinite nuclear masses 
Ml = M2 = 00, and obtain therefore a parametrical dependence of the Hamiltonian on R. 

The terms Hen, Hee, Hnn denote two-particle operators describing the classical Coulomb 
interaction between electrons (e) and nuclei (n). The term H°' in ([1]) is a one-particle 
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operator considering the kinetic energy of electrons of mass m, with A, referring to the 



Laplacian of the i-th electron. This form o: 
that attempts to include relativistic effects 



' the operator represents the simplest model 
S-sl- It does not give accurate numerical 



agreement with experimental observations, however, it does provide a qualitatively sensible 
description. For a = we retrieve the nonrelativistic limit of the kinetic energy of the i-th 
electron of — |A. The simple model chosen here has in common with all other relativistic 
approaches that it predicts instability of large atoms and molecules. The relevant parameter 
in this context is the product of the atomic numbers Z^,, k = 1,2, and the fine-structure 
constant a. The correct critical value of Zka as derived from the study of the Dirac 
equation is expected to assume Zka = 1. If the interest lies in the consideration of the 
limiting behaviour for ^ oo, one is therefore required to bound Zka by imposing 
a — )■ 0, as was done in 1^]. While this is from a physical perspective not sensible as the 
fine-structure constant exhibits an experimentally established value of approximately 1/137, 
it permits though a mathematical consideration of the asymptotics. The above discus- 
sion referred to the case of a diatomic system; the same is true for more complicated models. 



Concerning the abovementioned relativistic instability quantitatively, we recognize that 
our simple relativistic model exhibits a known deficiency in that its point of instability 
is in contradiction with the expected value of Zka = 1. In our model, the ground 
state energy E{N,Z,R:a) is finite if max {Zka} < 2/7r, but E{N, Z, R; a) = — oo if 
max {Zka} > 2/7r (see pHsl, as well as y for further discussion). This imposes that 
the atomic number has to be smaller than or equal to 2/7ra ^ 87, thereby contradicting the 
observation of atoms with atomic numbers larger than 87 to be stable. It is therefore clear 
that while our model is qualitatively reliable, we cannot expect good quantitative agreement. 



AT 



The operator H{N,Z,R;a) acts on the fermionic space /\L^(R^;C^). The energy 

i=l 

E{N, Z, R; a) as a function of R is given as 

E {N, Z, R; a) = inf spec H{N, Z, R, a). 
We introduce the energy minimum (or rather infimum) of the molecule 



E{N,Z;a) = ini E{N, Z,R;a). 

R 
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That E{N,Z;a) > — oo, i.e., the stabihty of relativistic molecules, was first proved by 
Daubechies and Liebj^. A more general proof is due to Conlonj^] and Fefferman and 
de la Lave[8|. 

We say that the molecule has a stable Born-Oppenheimer ground state if the following two 
requirements are satisfied. 

(1) - 

E{N, Z; a) < E{N, Z,R = 00; a) := lim E{N, Z, R; a). (2) 

(2) . The infimum is attained for some Rqi i.e., 

E{N,Z;a) = E{N,Z,Ro;a), (3) 

with the further requirement that it is an eigenvalue below the essential spectrum for 
H{N, Z,R;a) if R = Rq. 

The first requirement ensures that the individual atoms stay bounded, the second require- 
ment ensures that all electrons remain bounded to the molecule when the bond distance is 
Rq. We identify in the physical context Rq as an equilibrium bond distance, i.e., the distance 
between Ri, R2 in a stable molecule. Note that we do not exclude the existence of more 
than one Rq; this case plays no role in the present consideration though. 

If stability holds, by which we mean the existence of a stable Born-Oppenheimer ground 
state as defined above, we introduce the binding energy for a molecule as 

AEb{N, Z; a) = E{N, Z,R = 00; a)- E{N, Z, Rq; a) > 0. (4) 

We state first the main results of this paper before establishing the proofs in the subsequent 
sections. The final proofs of Theorem 1 and Theorem 2 can be found in Sec. V and Sec. VI 
respectively while Sees. II-IV introduce necessary prerequisites. 

Theorem 1 Assume H{N, Z, R; a) has a stable ground state for some Z = Zi + Z2, with 
max{Zka} < 2/7r. Let e > be such that max{Zka} < -{1 — e). Then, independent of 

k fc 

particle symmetry. 



Z1 + Z2- \ 2 V 4 



<N l + Jl + Sa{e,r)], (5) 



with a{e, t) = 2[1 + e ^ {1 — e) ^]t, where t > is the constant in Lemma 3, lf2J\) . 



Assuming fermions, we obtain instead 



Z1 + Z2- I 2 V 4 Ar2/3 



'''' .Wi + .i + ^|. (6) 



A consideration of the case Zi ^ Z2 shows that the bound on N is actually controlled by 
the smaller of the two atoms. 



An upper bound on is given by Lieb 10|] as A^ < 2(Zi + Z2) + 2, valid for the nonrelativistic 



as well as the relativistic case. In the latter, Dall'Acqua, S0rensen and Stockmeyer 111 ] 
proved this bound to hold for max {Zka} < 2/7T. 
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Theorem 2 There exist constants ci, C2 > such that if H{N, Z, R; a) has a stable ground 
state on the fermionic space for Z = Zi + Z2, Z = N, with max{Zka} < 2/7r, A; = 1, 2, then 

k 

0<AEt{Z;a)<CiZ^-^/''^ (7) 
i?o>C2^-^/'+"/'^°. (8) 



From Thomas-Fermi theory we know that the bulk of electrons around one nucleus is at a 
distance Z~^^^ from that nucleus. Thus the theorem states that the internuclear distance 
is much larger than the radius of the bulk electron cloud. 



The approach in the proof of Theorem 1 consists of the introduction of localization functions 



xi,X2 e 



satisfying 



Xi{x, RY + X2{x,Ry = l. 



(9) 



The nonrelativistic terms ifee + H^n + Hnn can be treated according to a localization 
argument following The nonlocal operator H°' will be treated separately, applying 
here results from 2|. We will begin with this consideration in Sec. [TTl before introducing 



the localization argument in Sec. Illli In Sec. |IV] we establish an estimate on Rq as 
a function of A^. We introduce in this context Lemma 3 and Theorem 4, the proofs 
of which can be found in the appendix. These results will be used in Sec. |V] to prove 
Theorem 1, i.e., to derive a global bound on A^. Theorem 2 will be proved in Sec. IVIl where 

n 

we will apply the recently proved relativistic Scott correction[2[ to derive bounds on AEb, Rq. 

We will in the following denote the wave function of the many-particle system by with 
* e A^^^(K^ X and = 1. 

II. RELATIVISTIC IMS FORMULA 

We will first consider the kinetic energy operator H°' = ^ separately by application 

i 

of the relativistic Ismagilov-Morgan-Simon (IMS) formula. This formula was proved in 
[2! generally for a family {C,u)u&m of positive bounded C^-functions on with bounded 
derivatives which, given a positive measure dfi on Ai, fulfill j\^C,u{xydf^{u) = 1 for all 
X G R^. We consider in the present context the set of two localization functions Xi, X2- The 
formula reduces in this case, for any ^ G H^^^{R^), to 

= (V^|XiT-xi|V^)-(^l^xJV^) 

+ {i^\X2T''X2\i^)-m^M, (10) 

with L^j, denoting an operator with integral kernel 

L^,ix,y) = {2n)-^a-' \x - y\'^ [a'^ \x - y\) [xfc(x) - xAvt , (H) 
where K2 is a modified Bessel function defined by 

/•oo 

7^2 (t) = t \ e-'-^hHs. (12) 
Jo 

Note that our choice of unit system differs from {2!, Theorem 13 {Relativistic IMS formula), 
in setting m = 1. 
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Employing the relativistic IMS formula, we can reformulate the following matrix element as 

i 

-mL^,ix,,ym - {^\L^,ix,,ym}- (13) 
III. LOCALIZATION ARGUMENT 

We consider now two-cluster decompositions 6 = {61,62) of {1, . . . , A^}. The interactions 
within the two separate clusters are described by Hg-^ , Hg^ , with Hg = Hg^ + Hg^ and 

while the intercluster potential Ig can be given as 



J 6 62 



and we recognize that H = Hg + Ig. 



The space of localization functions can be correspondingly regrouped as 

N 

n ixi{x.f + X2{x,f) = 5^ n n ^^i^.o^ = e 

i=i 9 ieei jee2 e 

with Xe = n Xii^i) n X2(a;j), omitting here and from now on out of reasons of simplicity 

iG6»i j£e2 

the dependence in R. 



We consider, starting from (IT^ . that 

TV 

(vl/|i7"|x[/) = { {^\Xi{x^) Xiixi) n ixii^j)' + X2(x,)') 1^) 

N 

+ {^\X2{x^) Tr X2{x^) n ixiixj)' + X2(a;,)') 1^)} 

i 

e i 



and furthermore reformulate the Coulomb terms straightforwardly 

e 
e 

Additionally, we recall that if H has a stable ground state then 

E{N, Z- a) = inf spec H{N, Z, R, «)|^^^^ < Ee{N, Z, Rq; a) 
= inf spec He{N, Z, R; 0!)\j^^j^^. 



This is the Hunziker-van Winter-Zhislin (HVZ) theorem 12j, in its formulation for this 
operator. 

We can therefore estimate 

R{^\H\^) = Rj2{{'^\XeH''Xem + ^Xe {H - xe^' 
e 

-RJ2{{'^\Lx.i^'.y^m + {mx.i^^,y^m\ 

i 

= RY,{{MxeHeXem + nrnhxlm 
e 

-RY,{mLA^^,y^m + mLxM,yrm\ 

i 

>RoY,{Eemxim + mhxim] 

e 

-R,Y,{('^\LAx^,ym + mLA^^,y^m] 

i 

>Ro M Ee + Roy^iWexll"^) 

e 
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We obtain a bound on the localization error {ilj\L^^{xi,yi)\ilj) in the following as 

J i-^ y\ 

<(27r)-2a-3 II Vxk \\l J mx)\my)\K2{a~'\x-y\)dxdy 
<(27r)-2a-3 II Vxfc IlL j \Hx)\^K2{a-'\x-y\)dxdy 



= (271)-' II Vxk \\L J m^rdx J K2{\y\)dy 

= I II Vx. IIL II^IIL (17) 

where in the third step we used the Cauchy-Schwarz inequality and K2 > 0. We obtain from 
this the final estimate 

-AE,{N,Z;a)Ro > -^iVi?o {II Vxi ||L + II ||L} 

N 



\xi-Ri\ \xi-R2\ 



N 



+ ^1 ^2- (18) 



IV. ESTIMATE ON i?o 

Following [l| we introduce explicit expressions for the localization functions which exhibit 
a dependence on the parameter fi = Z2/Z1. We impose for convenience Zi > Z2, therefore 
/i < 1. Introducing 

^ = ^ + 7^7 — r4T(^i - -R2) and R= — |— -(i?i-i?2) 
2(/i + l) n + l 

we obtain 

x + nR = x — R2 and a; — i? = a; — iii, 

using that i?2 = — -Ri- 
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We define the localizing functions as 

\x + fiR\ 



^\x-R\ 

X2{X) = _2 , , 

,/jrTi{\x\^ + fiRy/^ 

with R = \R\. It is straightforward to see that 

Xii^y + X2ixy = 1. 

We evaluate 



< ^J^R-\ (19) 
A* 

We will use this expression for R being the equilibrium bond distance Rq. We recognize that 
an estimate on Rq^ is necessary to subsequently deduce an estimate on the localization error. 
Obtaining this estimate is more involved in the relativistic case than in the nonrelativistic 
one due to the presence of the critical point if Z^a = ^ for some k. We start out by 
comparing the energy E of the diatomic system with the situation of localizing all electrons 
on one of the centers, say Ri. We obtain using the stability conditions (|2]), (|3]) that for the 
energy i?"* of a single atom 

E''\N,Zi;a) > lim E{N,Z,R;a) 
>E{N, Z,Ro;a) 

>E^\N,Z^ + Z2,a) + ^. (20) 

The last inequality, the lower bound on the energy of a united atom, is simple to derive for 
the nonrelativistic case, as was done in 



13| (Sec. (4.6.14)). The extension to the present 
case of the relativistic operator is straightforward as the form of the one-particle operator 
is identical for the molecular case and the united atom. 



Establishing an upper bound on E"'^{N, Z\\ a) — E°-*(N, Z\ + Z2; a) will allow us to deduce 
an upper bound on -Rq ^• 
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Since Z ^ i?"*(A^, Z] a) is a nonincreasing, concave function we have 



E-\N, Zi; a) - E^\N , Z^ + Z^; a) 
where [^■]_ refers to the left derivate. 
We note that if 

dE''\N,Z;a) 

we trivially have a lower bound of zero while if 

'dE"'\N,Z;a) 



> 



dE^\N, Zi + Z2;a) 
dZ 



dZ 



<0, 



then E''\N,Z;a) < 0. In this case there is an 1 < n < such that E''\N,Z;a) 

E"'^{n, Z] a) < E^*{n — 1, Z]a) and there is an n-particle wave function such that 

H-\n,Z;a)^n = E-\n,Z;a)^n, 

with H°-*{n, Z; a) being the Hamiltonian of the atomic system. 



By the Feynman-Hellman Theorem we have 



dE''\N,Z-a) 
dZ 



and have therefore reformulated the problem into establishing an upper bound on 
ELi R^n) when Z- a) < 0. 



This we will achieve in the following by first providing a lower bound on 



E 

i=l 



E 



1=1 



rpa 



for Za < I (Lemma 3, fl2T]) ) which will be accomplished by use of the combined Daubechies- 
Lieb-Yau (DLY) inequality (see furthermore js], Q for Daubechies inequality and 
Lieb-Yau inequality). This lower bound will then be used in the proof of Theorem 4, 
namely the upper bound on ^\l/„, X^ILi jF| ' '^^^^ — f (-'- ^ ^ft^r having 

established Theorem 4, we are in the position to return to ( EUj) and provide a bound on 
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Rq \ We recognize that our approach fails for the critical case of Za = -. 



Lemma 3 If Za < - then for all n > 1 

— r — 

z 



E 

i=l 



Xi 



> -Z^T 



n without particle symmetry 
n^^^ assuming fermions 



-Z^M 



(21) 



for some constant r > 0, with the function K{n) defined by lf2T]) . 



Theorem 4 Assume Za < ^{1 — e), < e < 1. Assume eigenfunction of H°'^{n, Z; a) 
with eigenvalue E"-\n, Z; a), where 1 < n < N and E°-^{N, Z\ a) = E'^\n, Z; a) < E''\n - 
l,Z]a). Then 

'dE^\N,Z-a) 



dZ 



> - f ^n, j > - [1 + (1 - ^y^] <ri) Z. (22) 



The proofs of Lemma 3 and Theorem 4 can be found in the appendix. 

We are now in the position to obtain our final estimate on Rq^ . Starting with the bound 
on the united atom we reformulate 



< E^\N, Zi; a) - E''\N, Z^ + Z^] a) 

< -Zo 



Rq 

'dE^\N, Z, + Z2;a 



dZ 

< Z2 {Zi + Z2) [1 + e-^ (1 - e)-'] k{N) 
<2 [1 + £-1(1 -£)-!] k{N)Z^Z2, 



with Z2 < Zi. We have used that n < N. 



We therefore obtain 



, , without particle symmetry ^ ^ 

Ro'<^ie,r) { . \ . (23) 

i\ ' assummg fermions, 

where a{e, r) = 2 [1 + e^^ (1 — e)^^] r, for some constants 1 > £ > 0, r > 0, with r intro- 
duced in Lemma 3, (pTl) . 
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V. GLOBAL BOUND ON 



The estimate on Rq as a function of allows us to estimate the localization error from f ll7p 

as 

RoML^,{x,ym\ + Ro mL^,{x,ym\ <^R^ {\\ Vxi \L + II Vx2 ||L} 

<3Ro II (Vxi)' + (VX2)' Hoc 
<3aiE,r)N ^^^'^^\ (24) 

Inspecting the remaining terms in f[T5]) . we recognize the difficulty in obtaining an estimate 
on the repulsion of electrons between different clusters and will neglect the term at this 
point. We thereby obtain from (fT8|) 

> -3aie,r)^J^N^ - i^o V(v&| ^^"^^ + + ^^^'T ^' 1^) + Z,Z,. (25) 



We estimate following 



{Z2\x + fiR\ + fiZi\x - R\Y < {Zl + Zlii) {\x + iiR\^+ \x-R\^) (26) 



= {Z^ + Ziii){ii + l)(|x|' + iiR 

and obtain with ^ = Z2/Z1 



2 Ar -2'i^2 , / ^1-2^2 ^ 



0>-3a(e,r)iV^-iV-^ + 

therefore 



Zi + Z2 - \^2 V 4 
We compare this result to the nonrelativistic case where Solovej[l| proved 

/^<iV(3/2). 

The difference in the bound on is attributable to the prefactor in the first term of f l25|l 
which is in the nonrelativistic case — | instead of — 3(j(£:,r). Aside of the more demanding 
assessment of -Rq^, which in the nonrelativistic case yields an estimate of Rq^ < | A^, we 
recognize two further steps in the present derivation which infiuence our bound. The simple 
estimate of || Vxi ||^ + || Vx2 ||L < 2 || (Vxi)^ + (Vx2)^ lU in the second step of 
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introduces a factor of two. A further factor of three has to be attributed to the estimates in 
ffT7|) . In particular, we bound the expression \xk{x) — Xk{y)\ by 

\Xkix) - Xkiy)\ <|| Vxfc Hoc -{x - y\ 

before employing subsequently the Cauchy-Schwarz inequality. Instead, it would be 
preferable to construct explicit forms for Xij X2 which would allow a tighter estimate 
of (tplLy^Jip) while at the same time permitting an advantageous evaluation of the 
electron-nuclear attraction term in f l26|) . With respect to the localization error, one 
has in particular to construct a form of localization functions which results in a can- 
cellation of the denominator \x — y\ in order to obtain a bound of | j-?/') | < 00. 
In the present context, we have decided to employ the particular forms of Xiy X2 given 
in l| which result in an algebraically simple expression of the bound on and recog- 
nize the scope that exists in the possible improvement on the bound of the localization error. 

Considering alternatively the fermionic bound on Rq^ in our estimate, we obtain 



Z1 + Z2 - I 2 V 4 A^2/3 

This can be compared to Benguria, Siedentop and Stockmeyer[15|] who find a very similar 
expression for the case of homonuclear relativistic molecular ions. Their proof holds for 
Za < 1/2, where Z = Zi + Z2 = 2Zi. 



VI. BOUND ON AEb, Rq 

The first rigorous investigaton of the limit Zk ^ 00 with Zk<y bounded was given by 



S0rensen[16|. The leading asymptotics of the ground state were established to be deter- 
mined by Thomas-Fermi theory. The first correction term, the Scott correction, and its 
dependence on Z^a, was proved by Solovej, S0rensen and Spitzer|3] for a neutral molec- 
ular system, i.e., N = Z = Zk, in the framework of the Born-Oppenheimer formulation. 



We state their main theorem, the relativistic Scott correction, for the case of a diatomic 
system. Let Z = Zi + Z2, z = (2:1,-22) = Z^^Z with zi,Z2 > 0, thereby zi + Z2 = 1. We 
define r = Z^^^R, with \r\ > tq for some tq > 0. Note that r = ri — r2, with ri, r2 E M^, 
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furthermore r = jri — r2|. Then there exist a value E'^^{z,r) and a universal (independent 
of z and r) continuous, nonincreasing function S : [0, 2/7r] — )■ R with 5(0) = 1/4 such that 
as Z — )■ oo and a — )■ 0, with max {Z^a} < 2/7r, we have 

E{Z, R; a) = Z'l'' E^^{z, r) + 2 Zl S{Z^a) + 2 ^(Zaa) + ^(Z^-i/^O). (27) 

The error term means that |(9(Z^~^/^'') | < cZ^"^/^*^, where the constant c only de- 
pends on tq. Moreover, the Thomas-Fermi energy satisfies the scaling relation 
E^^{Z, R) = Z^/^E^^{z, r). We call E^^{Z, R) the Thomas-Fermi energy of the molecule. 



We consider from 



that 



Ro > 



.iV-l/3 



cr{e,T 



z-^i\ 



therefore asserting that Rq > Z ^^^tq for some constant ro > 0. 

We will use in the following evaluation scaling properties from Thomas-Fermi theory. In 
particular, it is established (see 17|] for a review of TF-theory) that 



E^^(Z,i?) = Z^/3^^^(z,r). 



(28) 



We introduce furthermore the estimate on the split o 
theory which was derived in {l| based on Brezis-Lieb 



' a diatomic system in Thomas-Fermi 

3 



E^^(z,r) > E^^iz,) + E^^iz,) + cr-\ 



(29) 



with some constant c. E'^^{zn) denotes the scaled Thomas-Fermi energy of a nucleus of 
charge n at the origin. 



We consider now the case of a diatomic molecule and obtain by application of the relativistic 
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Scott correction (127|) and furthermore (1251) . (12^ that 

H{N,Z,R;a) > E{Z,R;a) 

> Z'/'E^'^izi) + 2 Zl S[Z^a) + Z^/^^'^f^^^^ ^ ^ ^2 ^^^^^^ 

= E^^{Zz^) + 2 + E^^^Zz^) + 2 Z| S^Z^a) 

> ^(Zi; a) + E(Z2; a) + ci r"^ Z^/^ - cq Z^-^/^o, (30) 
with constants cq, Ci (note that the definition of cq changes in the last step). 

We conclude 

E{N, Z, i?; a) - E{Zv, a) - ^(Zs; a) > ci i?o ^ - co Z^'^/^", (31) 
and can therefore estimate 

0>cii?o'-coZ2-i/30 
i?o>cZ-i/-^+ii/2io (32) 

The constant c is independent of Z and i? and may assume different values in different 
inequahties. 

We compare to the nonrelativistic case where Solovej [ll proved 
for < Z and e = 1 /TO. 

We note that the reason for our bound to be better than the bound of Solovej is that we 
used that the energy is known up to the Scott correction. The same can be done in the 
nonrelativistic case and would give a similar bound. The nonrelativistic Scott correction 
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was proved by Ivrii and Sigaljl^ after the publication of 
Benguria, Siedentop and Stockmeyer 



151 ] state for homonuclear relativistic molecular ions a 
bound on Rq (for Za < 1/2) which, translated into our unit system, corresponds to an error 
of the order Z~^. The difference to our estimate reflects their consideration of boltzonic 
electrons. 

Recalling that in the case of a stable ground state 

AEbiN, Z; a) < E{Zi; a) + E{Z2] a) - E{N, Z, R; a), 
we deduce furthermore a bound on AEb{N, Z; a) of 

AEh{N,Z;a) <cZ^-'/^^. (33) 

VII. ACKNOWLEDGEMENT 

I would like to thank Jan Philip Solovej for providing the framework for this work as well 
as for numerous insightful and stimulating discussions. 

This work was supported by a fellowship within the Postdoc-Program of the German Aca- 
demic Exchange Service (DAAD). 

Appendix A: Proofs 
1. Proof of Lemma 3 

We begin with a separation of the nucleus-electron attraction by considering a ball |x| < r 
for some r to be chosen later. The contribution from the outside region is estimated by — ^. 
We assume n > 1. 



i=l 



i=l 



Xi 



— , Al 

r 



where l{|2;|<j.} is the characteristic function of the ball < r}. 
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A lower bound on the right-hand side is in the fermionic case obtained through a sum over 
all negative eigenvalues. In the case of no particle symmetry, a lower bound is obtained by 
approximating the bosonic states through n times the lowest eigenvalue, which is then in 
turn estimated by the sum of all negative eigenvalues. 



z 




nZ 




l{|x.|<r} 


r 


z 




nZ 




l{|x,|<r} 


r 



> 2Tr 



z 




nZ 


\x\ 


l{|xi<r} 


J. 


z 




nZ 


\x 


l{|x|<r} 


J, 



for spin 1/2 fermions(A2) 
without symmetry. (A3) 



E 

i=l 
n 

T^-^U.A<rri - — >nTr - 

i=l 

We will use the combined Daubechies-Lieb-Yau inequality {2!, [l^ to finalize the proof of 
Lemma 3. Note again that our form differs from [2], Theorem 16 [Combined Daubechies- 
Lieb-Yau inequality), in setting m = 1. 

We state the combined DLY inequality in its application to the simpler case of a single atom. 



Assume a function W G Lj^^CR^) which satisfies 



W(x) > —- — r — Cva when \x\ < a, 

\x\ 



with au < 2/7r, a > 0, and a constant C. Then 



Tr [T° + > -Cz/^/'a^/' - C J \W{x)^f/^d 

a<\x\ 

-Ca^ j \W{x)_\^dx, 



a<\x\ 



where, as above, T" = a/— a ^ A + a ^ — a and = —A/2 when a = 0. 



(A4) 



(A5) 



With z/ = Z, we can evaluate 
Tr 



7 

- ^1 

\x\ 



{\x\<r} 



> -CZ'/'a'/'-C J (^)'^'^^ 

|a;|<r 

-Ca^ / (^) 

a<|a;|<r 

= -CiZay/^Z^ - CZ^/V^/^ _ c(ZafZ^ 



(A6) 



18 



where we have used that Za < -. 



With the choice of r of 



r — Z ^r?!"^ for fermions (A7) 
r = Z"-*- without symmetry, (A8) 



and using that n > 1, we complete the proof of Lemma 3. 



2. Proof of Theorem 4 

Using Lemma 3 we can state 

1 • 1 r^* 

We obtain an estimate on 'YTi=\ T^^n) through an assessment of the energy of a single 
atom 

0> = (AlO) 

n n n „ 



> £(^„,^7;"^„)-(l-£)«;(n)-— (All) 

1=1 ^ ' 

using again Lemma 3 in the last step, asserting that -^ot < | by the assumption of Theo- 
rem 4. We can conclude that 

Y T^'^n) < ^ ^/^N Z\ (A12) 



=1 



and thereby complete the proof of Theorem 4. 
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